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1.  INTRODUCTION 

Financial  engineering  relies  on  an  approximation  —  that  the  value  of  stocks,  bonds,  curren¬ 
cies,  etc.,  follow  either  a  continuous  path  or  a  binary  tree.  The  methodology,  going  back  to  Black 
and  Scholes  (1973),  Harrison  and  Kreps  (1979)  and  Harrison  and  Pliska  (1981)  —  see,  e.g.,  Duffie 
(1996)  for  a  contemporary  overview  —  has  been  a  powerful  factor  in  the  development  of  options 
markets.  And  yet,  the  cracks  in  the  theory  show  up  in  phenomena  such  as  the  volatility  smile  (see, 
for  example,  Ch.  8.H  (p.  180-185)  of  Duffie  (1996)). 

Most  attempts  at  tackling  the  deviations  between  theory  and  reality  have  rehed  on  creat¬ 
ing  additional  flexibility  in  the  continuous  path  model,  in  particular  through  stochastic  volatility 
models  (including  HuU  and  White  (1987),  BaU  and  Roma  (1994)  and  Duan  (1995))  and  implied 
volatility  trees  (Dupire  (1994),  Derman  and  Kari  (1994),  Rubinstein  (1995)).  Another  Une  of 
research  has  abandoned  continuity  (or  trees)  altogether,  and  instead  focused  on  equilibrium  con¬ 
siderations  and  on  approximate  hedging  (Merton  (1976),  FoUmer  and  Sodermann  (1986),  Aase 
(1988),  FoUmer  and  Schweizer  (1990),  Schweizer  (1990,  1991,  1993,  1994),  EUiott  and  FoUmer 
(1991),  ColweU  and  EUiot  (1993),  Mercurio  and  Runggaldier  (1993),  to  mention  some). 

There  are  advantages  to  both  these  approaches.  Accepting  discontinuity  leads  to  a  more 
accurate  model  of  prices,  which  typically  move  in  multiples  of  fixed  increments.  It  also  permits  a 
more  useful  description  of  shocks  such  as  stock  market  crashes,  currency  devaluations,  and  similar 
occurrences.  On  the  other  hand,  discontinuity  carries  with  it  an  abandonment  of  the  arbitrage 
based  theory  of  pricing.  Even  in  a  fiexible  continuous  path  model,  such  pricing  can  be  used,  if 
need  be  by  including  coefficients  representing  the  market  price  of  suitable  untraded  securities.  The 
advantage  of  continuity  is  pragmatic,  but  no  less  compeUing  for  that  reason. 

The  purpose  of  this  paper  is  to  present  an  intermediate  strategy,  which  can  draw  on  the 
strength  of  both  the  approaches  described  above.  We  shall  suppose  that  the  stock  price  moves  in 
jumps  of  \arying  size,  and  then  let  the  intensity  of  the  jumps  increase  as  their  size  decreases.  The 
stock  price  wiU  then  converge  to  a  geometric  Brownian  motion,  but  it  wiU  also  be  possible  to  write 
an  asymptotic  expansion  for  the  price: 


price  =  Black-Scholes  price  -|-  -f-  o(^^/^). 


(1.1) 
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Here,  ^  is  a  quantity  which  measures  the  distance  between  the  actual  stock  price  model  and 
geometric  Brownian  motion. 

In  this  setup,  .R  is  the  price  of  a  residual  security  in  a  limiting  security  market,  and  it  involves 
observable  constants  and  the  market  prices  of  a  low  dimensional  set  of  untraded  securities.  To  first 
order,  therefore,  captures  the  cost  of  discontinuity,  but  at  the  same  time,  it  is  obtainable 

from  arbitrage  considerations. 

It  would,  obviously,  be  more  desirable  to  have  the  true  price  rather  than  an  asymptotic 
expansion  for  it.  (1.1),  however,  can  be  seen  as  an  analysis  of  the  sensitivity  of  the  Black-Scholes 
price  to  small  deviations  from  the  standard  model,  and  as  such  it  should  be  a  useful  tool  both 
theoretically  and  practically.  -  Note  that  compared  to  the  results  in  WiUinger  and  Taqqu  (1991). 
DufRe  and  Protter  (1992),  Cutland,  Kopp  and  WiUinger  (1993),  and  Amin  and  Khanna  (1994). 
(1.1)  provides  not  only  a  limiting  result,  but  also  the  first  order  remainder  term.  In  this  sense,  we 
are  closer  to  Foster  and  Nelson  (1996)  and  similar  work  by  both  these  authors.  The  remainder 
term  is  critical  in  the  further  analysis. 

In  addition  to  providing  a  more  accurate  modeling  of  stocks,  discontinuity  also  permits  a 
realistic  description  of  statistical  uncertainty.  This  feature  is  not  reaUy  available  in  the  continuous 
time  model.  Under  continuity,  the  volatility  is  observed.  Hence,  if  it  is  constant  there  can  be  no 
statistical  uncertainty  (since  the  drift  does  not  have  an  impact  on  the  price).  One  can  go  through 
various  attempts  to  induce  such  uncertainty  by  making  the  volatility  nonconstant,  but  this  wiU 
only  recreate  the  problem  at  a  meta  level.  Consider,  for  example,  the  case  of  endogenous  variation. 
It  is  highly  desirable  to  model  this  variation  through  a  stochastic  volatifity  model,  say, 

=  u{a^)dt  +  ^{al)dWt,  (L2) 

The  reason  for  this  is  that  it  is  difficult  to  model  mean  reversion  with  a  non  random,  time  varying 
(7^ ,  and  mean  reversion  is  widely  believed  to  be  a  feature  of  the  volatility  process.  The  only  possible 
statistical  uncertainty  would  then  reside  in  u  and  7.  The  latter,  however,  is  observed  on  the  domain 
where  cr^  has  actually  varied  (and  is  completely  unknown  otherwise),  and  u  as  such  does  not  affect 
the  price  of  options,  as  it  has  to  be  modified  by  the  market  price  of  risk  in  Also,  v  is  practicaJlv 
observable  in  the  long  run  (if  one  incorporates  information  on  options  prices;  see  Renault  and 
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Touzi  (1996)).  One  can  make  v  and  7  time  varying,  but  then  the  question  of  incorporating  mean 
reversion  comes  up  again.  And  so  on. 

One  can  also  model  as  a  function  of  5,  or  in  various  other  ways,  but  as  far  as  the  author 
is  aware,  there  is  no  known  way  of  escaping  the  problem  we  have  described  in  the  framework  of 
the  continuous  path  model.  One  can  get  around  it  by  assuming  that  the  volatility  is  observed  by 
market  participants  but  not  by  econometricians  (Renault  (1995)),  but  this  is  not  a  particularly 
satisfying  solution.  One  can  also  adopt  the  approach  of  Gottlieb  and  Kalay  (1985),  Ball  (1988),  and 
Kleive  (1993)  of  thinking  of  stock  prices  as  rounded  or  truncated  values  of  a  continuous  process, 
but  again  this  raises  the  question  of  whether  stocks  can  be  hedged  continuously  or  not. 

This  is  not  to  say  that  the  problem  of  inference  in  models  for  the  behavior  of  financial  instru¬ 
ments  is  not  a  difficult  one.  A  substantial  literarture  exists,  see,  e.j.,  BoUerslev,  Chou  and  Kroner 
(1992)  for  an  overview  of  ARCH  and  GARCH  models,  and  Dachuna-CasteUe  and  Florens-Zmirou 
(1986),  Bibby  and  Sorensen  (1994,  1995),  Danielsson  (1994),  Hansen  and  Scheinkman  (1995),  Ait- 
Sahalia  (1996),  Kessler  and  Sorensen  (1995),  Renault  (1995),  Foster  and  Nelson  (1996),  Gallant 
(1996),  Hansen,  Scheinkman,  and  Touzi  (1996),  Pastorello  (1996)  and  Renault  and  Touzi  (1996) 
for  a  selection  of  important  recent  contributions  to  inference  in  stochastic  differential  equations 
(SDEs)).  It  is  just  that  there  is  little  integration  between  the  theory  of  inference  and  the  theory 
of  pricing. 

Arguably,  statistical  uncertainty  is  not  a  problem  —  one  can  just  use  implied  quantities. 
Things  are,  however,  not  quite  so  simple  in  the  presence  of  many  unknowns,  as  we  shall  see  in 
Section  3.  First,  however,  we  describe  the  model  and  derive  the  expansion. 

2.  ASYMPTOTIC  ANALYSIS 

2.1.  The  Model 

We  suppose  that  log5t  is  a  pure  jump  process  whose  jumps  have  intensity  Aj,  and  that 
the  jump  at  time  t  has  (random)  size  This  can  formally  be  carried  out  by  letting  logS*  be 
the  integral  of  the  Zt  process  with  respect  to  the  random  measure  that  has  intensity  At  (see,  e.g., 
Ch.  Il.l.d  (p.  71-74)  of  Jacod  and  Shiryaev  (1987)).  The  asymptotics  are  then  implemented  by 
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letting  the  intensity  go  to  infinity  while  the  jump  size  goes  to  zero.  To  be  precise,  let  St  =  be 
indexed  by  a  <5  which  tends  to  zero,  and  suppose  that  the  first  four  cumulant  variations  of  log  5/ 
(see  Section  6.1  (p.  30-32)  of  Mykland  (1994))  exist  and  have  the  form  Jq  fisds  =  fiit  +  Op(l), 
Jo  ^^ds^  +  Op(<5^/^)  and  +  Op((5).  We  also  assume  that  £’sup(Alog5t)^  =  o(^),  which 

in  particular  implies  that  higher  order  optional  variations  of  log5t  are  of  order  Op(<5). 

The  conditions  assure  that  log5‘^  converges  to  a  Brownian  motion  with  drift,  and  the  price 
of  an  option  wiU  be  shown  to  be  on  the  form  (1.1). 

Before  finding  the  price,  however,  we  shall  investigate  how  to  model 

2,2.  Volatility  Structure 

It  wiU  be  shown  is  Section  2.3  that  the  quantity  R  in  (1.1)  is  the  price  of  a  derivative  security 
in  the  limiting  probability  measure  obtained  when  6  goes  to  zero.  In  order  to  obtain  a  meaningful 
expression  for  ii,  this  limiting  security  market  must  have  an  information  structure  which  is  sensible 
for  our  purposes.  This  imposes  requirements  on  the  volatility  structure. 

First  of  all,  cannot  follow  (1.2),  since  that  would  lead  it  to  be  observable  in  the  limit.  In 
fact,  if  cr£  is  the  long  run  average  volatility,  one  needs  =  Op((5^/^)  to  avoid  this  problem, 

because 

[log5,log5]t +  (2.1) 

Jo 

50  otherwise,  <7^  would  be  observable  in  the  limit.  We  shall  therefore  suppose  that 

=  4  +  <■«,(-)  +  (2,2) 

For  similar  reasons,  we  let  A£,(5~^  be  the  long  run  average  intensity  process.  It  is  natural  to  suppose 
that 

At  =  (5"U£,  +  (2.3) 

(One  could  also  work  with  the  case  where  A(/^  tends  to  zero  or  infinity,  but  this  is  less  interesting 
as  one  factor  will  then  dominate). 

To  accomplish  (2.2)-(2.3),  one  can  suppose  that 

d{c'i,6\ty  =  u{a'l,6\t)dt  +  6^l^^{a1,6\t)dWt. 


(2.4) 
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Vt  =  then  takes  the  form 

dr)t  =  v'qtdt  +  'fdWt,  (2.5) 

where  i/'  =  i^'{crj^,XL)  and  7  =  7(<r£,  A^),  and  where  it  is  assumed  that  u{a\,\i)  =  0.  Since 
correlation  between  the  two  components  of  the  Brownian  motion  W  can  be  moved  into  the  7,  we 
shall  suppose  these  components  to  be  independent  standard  Brownian  motions.  Note  that  (rjt) 
is  mean  reverting  if  det(7)  0  (at  and  i/'  is  negative  definite  (in  the  sense  that  the 

eigenvalues  A  of  i/'  satisfy  Re(A)  <  0). 

The  principal  auxiliary  processes  in  our  system  are  then  the  observed  standardized  square 
variation 

=  d~^/2([log5,log5]<  -  (2.6) 

and  the  observed  standardized  activity  level 

nf  =  jumps  from  0  to  t  -  6~^XLt)  (2.7) 

These  quantities  are  ‘observed’  on  the  principle  that  to  first  approximation,  long  run  quantities 
like  a\  and  Xi  are  observable. 

Set 


One  then  gets  the  limiting  behavior  of  our  observables. 

Theorem  1.  Assume  the  conditions  of  Section  2,1.  Also  assume  (2.4),  that  7  is  continu¬ 
ous  and  that  its  eigenvalues  are  bounded  away  from  zero  and  infinity,  and  that  u  is  continuously 
differentiable.  Then  M  is  nonnegative  definite,  and  (2.2)~(2.3)  holds  with  respect  to  the  supre- 
mum  distance  and  with  rj  given  by  (2.5).  Furthermore,  ,W)  converge  jointly  in  law  to 

(S,^,  n,  IT),  where 

dSt  =  (/iL  +  Stdt  +  aiStdBt,  (2.9) 

6  = 

<^L 


(2.10) 
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and 


rj^du  + 


(2.11) 


where  /g  symmetric  square  root  ofM,  Vt  =  and  (5t,  VF/^*) 

is  a  standard  Brownian  motion.  The  convergence  in  law  is  in  the  space  of  cadlag  functions  on  [0,  T\ 
with  respect  to  the  supremum  distance. 


Proof.  See  the  Appendix. 


3.S.  The  Price  of  Options 

Suppose  we  wish  to  hedge  a  European  option  which  expires  at  time  T,  and  for  simplicity 
assume  that  the  risk  free  interest  interest  rate  r  is  constant.  K  C{S^t^a)  is  the  Black-Scholes  price, 
a  natural  hedging  strategy  is  to  use  a  “delta”  of  the  form  Set 

xf*  =  CiSo,0,cTL)+  f  Cs{St_,u,ai_)dSt\  (2.12) 

Jo 

where  =  e  .  If  we  let  Xf  =  then  by  numeraire  invariance,  there  is  a  self-financing 

strategy  for  X^,  the  price  of  which  is  C(5o,0,(ri). 

We  wish  to  evaluate  the  residual  exposure  to  risk  after  approximately  hedging  C{Sj,T)  with 
Xj..  Set 

4  =  <5-1/2  -  Xf)  ,  (2.13) 

where  C*  =  e  To  calculate  the  residual  price  R  in  (1.1),  we  find  the  limit  Rt  of  Rf,  and  from 
this, 

R  =  E%Rt\Xo),  (2.14) 

where  P  is  a  risk  neutral  probability  to  be  discussed  (cf.  Harrison  and  Kreps  (1979)). 

Before  we  proceed  any  further,  note  that  R  is  independent  of  hedging  strategy  of  since  the 
price  of  Xy  is  independent  of  this  choice.  For  the  purpose  of  evaluating  R,  therefore,  we  shall, 
without  loss  of  generality,  set 


df  =  OL. 


(2.15) 
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We  shall  use  C{S,t)  =  C{S,t.,ai).  This  does  not  mean,  however,  that  the  level  of  uncertainty 
about  at  does  not  have  an  impact  on  the  price  R,  as  we  shall  see  in  the  following.  First,  however, 
the  limit  result. 

Theorem  2.  Assume  the  conditions  of  Theorem  1,  and  also  that  CsssiS,t)  exists  and  is 
continuous,  uniformly  in  t,  for  5  €  A,  where  R  —  A  has  arbitrarily  small  Lebesgue  measure.  Then, 
(Rt)  converges  in  law,  jointly  with  and  in  the  same  topology  as  the  processes  in  Theorem  1,  to 

—  2  f  Ssi^ui‘f^')d^u 
Jo 

1  /■* 

+ fiSu,u)du,  (2.16) 

where 

f{S,  u)  =  3S^C*ssiS,  u)  +  S^ClssiS,  u).  (2.17) 


Note  that  the  conditions  of  Theorem  2  are  trivially  satisfied  for  call  and  put  options. 
Returning  now  to  the  price  of  the  residual,  note  that 

R  =  ^hT  js  (l  -  5o'Css(5o,0)  +  5'o"Csss(5'o,0)|  +  R'  (2.18) 

where 

R'  =  1^*  SlC^ssiSu,  u)dl  I  .Fo)  .  (2.19) 

This  is  because 

E-  (^f^C-CS,,!)  I  ;F„)  =  5;^C-(S„,0),  (2.20) 

and  because,  by  (2.10), 

=  + 

where  B*  is  a  Brownian  motion  under  P*. 

To  find  R',  we  need  the  value  process  of  In.  the  most  general  form,  suppose  that  Vt  and 
Wt  are  the  (two-dimensional)  market  prices  of  Vt  and  Wt,  respectively,  so  that  Vt  -  Vt  and  Wt  -  Wt 
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are  P* -martingales.  It  is  then  easy  to  see  that 


R'  =  \E*  SlC'ssiSu,  w)  -r  (1, du  |  Pc 


(2.22) 


where  rjt  is  the  solution  of 

drjt  =  i/'rjtdt  -I-  '•fWtdt 
Vo  =  Vo- 


(2.23) 


A  substantial  simplification  is  obtained  by  assuming  either  that  (ut)  and  (wt)  are  independent  of 
(St),  or  that  they  follow  a  linear  differential  equation 


d{yt,Wt)  —  AQ[t)dt  -f-  Ai(^t)(^vt,wt)*dt  -t-  A2{t)dZt,  (2.24) 

where  the  A’s  are  nonrandom  and  where  Zt  can  have  Bt  as  one  component,  with  the  other  com¬ 
ponents  independent  of  (Bt).  In  this  case,  it  is  straightforward  to  see  that 

R'  =  \S^Css(So,  0)E'(^t  I  Po).  (2.25) 


In  this  case,  since 

C;(S,  t,  a)  =  {T-  t)Css(S,  t,  a)S\  (2.26) 

one  can  alternatively  write  (1.1)  as 

price  =  6'(5'o,0,aft)  +  1^3^5^55(50, 0,crL)  +  o(S^^%  (2.27) 

where 

<^l  =  <rl  +  d'/'l  (^3  (l  -  +  IP*  (Jt  I  Po))  .  (2.28) 

Though  the  given  solution  features  skewness  (as  in  k^)  more  prominently  than  kurtosis,  the  latter 
is  present  through  k4  in  the  matrix  .1/  (in  (2.8)).  If  one  were  to  go  to  a  second  order  expansion, 
k^  would  also  be  present  in  other  ways. 
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3.  ESTIMATION 

3.1.  Implied  Quantities 

Three  cases  can  occur: 

i)  There  are  exchange  traded  options,  and  one  wishes  to  price  an  option  with  the  same  maturity 
as  an  exchange  traded  one. 

ii)  There  are  exchange  traded  options,  and  one  wishes  to  price  an  option  with  a  different  maturity 
than  the  exchange  traded  ones. 

iii)  There  are  no  exchange  traded  options  on  the  security  S. 

Reality,  obviously,  is  more  complex;  the  exchange  traded  options  could  be  American,  or  one 
might  wish  to  price  non  European  options.  For  simplicity,  however,  we  continue  to  concentrate  on 
European  options. 

Case  (i)  is  the  simplest  —  one  can  find  £'*(^r  I  -^^o)  by  regression  of  the  exchange  traded 
options.  Or  one  can  use  (2.27)  to  find  for  each  strike  price  after  correcting  for  the  second  term 
on  the  r.h.s.  of  (2.27).  If  the  approximation  is  good,  this  yields  a  fairly  flat  curve  for  ah  as  a 
function  of  the  strike  price.  One  should  not  expect  a  perfect  fit,  for  there  can  be  error  terms  of 
order  0(6).  How  to  carry  out  the  regression  or  the  averaging  of  the  cr/i’s  (scale?  weights?)  is  an 
interesting  question,  similar  to  issues  already  debated  in  connection  with  implied  volatilities.  We 
shall,  however,  not  go  into  this  issue  here. 

Case  (iii)  means  that  one  has  no  data  with  which  to  determine  the  relevant  market  prices 
of  risk.  In  this  case,  one’s  best  bet  would  probably  be  to  use  considerations  along  the  lines  of 
Schweizer  (1992)  and  Pham  and  Touzi  (1996). 

Case  (ii)  is  the  most  interesting  one  from  our  perspective.  One  has  observations  (or  estimates) 
of  E*{^j  I  JTq)  for  a  couple  of  expiration  dates  T,  and  one  wishes  to  interpolate  between  them. 
There  are  two  strategies: 

•  using  implied  quantities  only,  or 

•  using  a  combination  of  statistical  estimates  and  implied  quantities. 
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One  can  also  draw  on  equilibrium  considerations,  but  we  shall  not  consider  this  angle  here. 

Relying  only  on  implied  quantities  is  at  best  unwise,  and  often  impossible.  Even  in  the 
benign  case  where  the  A,’s  in  (2.24)  are  independent  of  t.  and  hence  approximately  observed  from 
long  run  data,  one  gets  that 

\  ^o)  =  fiVo^'>^o,^o,T)  (3.1) 

for  some  function  /.  Since  each  of  rjo,  vq  and  wq  are  two  dimensional  and  unobserved,  /  is  a 
function  of  three  to  six  unknowns,  depending  on  the  matrices  specifying  the  system. 

This  means  that  one  needs  three  to  six  different  option  expiration  dates  to  infer  the  quantities 
needed.  (Different  strike  prices  for  the  same  maturity  do  not  give  additional  information  for  the 
model).  Many  options  do  not  even  have  that  many  maturities.  Even  if  they  do,  this  is  a  dangerous 
reliance  on  the  validity  of  the  model.  It  would  be  safer  to  use  additional  observation  points  T  to 
assess  the  correctness  of  the  predictions  of  the  model. 

For  this  reason,  we  shall  in  the  following  taJce  a  statistical  approach  to  reducing  the  number 
of  implied  quantities  needed. 


S.£.  Uncertainty  Equilibrium 

Our  limiting  system  consist  of  observed  processes  5.  ^  and  n,  and  unobserved  processes  i] 
and  V .  It  seems  natural  to  assume  that  one’s  uncertainty  is  in  equilibrium,  in  the  sense  that  for 
all  t. 

Ft  =  r,  (3.2) 

where  Ft  is  the  posterior  covariance  matrix  for  given  the  data  up  to  t.  Alternatively,  one  could 
assume  that 

posterior  c.d.f.  of  at  t  =  F{x  -  ijt),  (3.3) 

where  F  is  fixed  and  where  is  the  (vector)  posterior.  The  reason  why  we  propose  this  is  that 
one  is  unlikely  to  systematically  gain  or  lose  information  about  the  unobserved  state  variables  over 
time.  In  other  words,  an  investor  will  typically  not  know  more  about  the  current  volatility  on 
March  22  than  on  January  7.  One  could,  obviously,  model  F  or  F  as  mean  reverting  processes. 
We  avoid  this  here  in  order  to  keep  the  model  as  simple  as  possible. 


11 


Uncertainty  equilibrium  has  implications  for  what  prior  and  posterior  distributions  must  look 

like. 

Theorem  3.  Assume  the  conditions  of  Theorem  1.  Also  suppose  that  v'  is  negative  definite 
and  that  M  is  positive  definite.  Then  the  equation 

TM-^T  -  u'r  -  r(i/')*  -  77*  =  0,  (3.4) 

has  a  unique  symmetric  positive  definite  solution.  Furthermore  the  following  statements  are  equiv- 
alent: 

(i)  (3,2),  and  the  distribution  on  tjq  given  the  data  at  time  0  is  Gaussian:  and 

(ii)  (3,3),  and  F  is  absolutely  continuous  with  full  support  on  B?, 

Under  either  of  these  equivalent  assumptions,  the  covariance  matrix  is  the  symmetric  positive  def¬ 
inite  solution  of  (3,4) •  Also,  the  conditional  expectations  satisfy 

dfit  =  u'fjtdt  +  TM-^{d{^t,nty  -  fitdt).  (3.5) 


Proof:  See  the  Appendix.  Note  that  a  reasonably  explicit  solution  to  (3.4)  is  constructed 
in  the  proof. 


In  other  words,  uncertainty  equilibrium  imposes  the  prior  variance,  and,  if  one  assumes  (3.6), 
the  entire  prior  distribution.  Under  these  circumstances,  Bayesian  inference  would  appear  to  be  a 
reasonable  choice  for  characterizing  the  available  information. 


In  the  above,  note  that  since  pi,  cr£  and  are  long  run  quantities  which  are  assumed  to  be 
observed,  and  Is  an  observed  process.  It  is  possible  to  take  but  we  then  suppose  the 

volatility  to  be  observed.  Otherwise,  the  assumption  of  uncertainty  equilibrium  would  be  violated. 
Now  set 


O'- /««■'“) 


(3.6) 


Since  f/i  is  a  conditional  expectation  given  the  data,  it  follows  that  Vf  is  a  martingale  in  the  filtration 
generated  by  the  observables.  Since  the  quadratic  variation  structure  of  {B,  V)  is  the  same  as  for 
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{B,V),  Levy’s  Theorem  (see,  e.g.,  Theorem  II.4.4  (p.  102)  of  Jacod  and  Shiryaev  (1987))  yields 
that  {B,V)  is  a  standard  Wiener  process  with  respect  to  the  observable  filtration. 

It  follows  from  Theorem  3  that 


dfit  =  u'fitdt  +  TM-'^^^dVt.  (3.7) 

Hence,  the  observed  system  is  characterized  by  (2.9)-(2.10),  (3.7),  and 

(6,nt)*=  frj^du  +  M^^^Vt.  (3.8) 

Jo 


3.3.  The  Effect  of  Uncertainty  on  Prices 

Formula  (2.19)  continues  to  hold,  but  we  now  have  a  new  observed  Doob-Meyer  decompo¬ 
sition  for  the  one  given  by  (3.8).  The  market  now  only  has  two  Brownian  motions  rather  than 
four,  and  if  one  takes  dVt  —  Vtdt  to  be  a  martingale  P*,  one  gets  that  (2.22)  remains  valid  with 


vt  =  Vt, 


(3.9) 


but  rjt  now  is  the  solution  of 

drjt  =  u'gtdt  -|-  TM~^^^Vtdt 
Vo  =  Vo. 

In  other  words  'ywt  gets  replaced  by  and  go  by  fjo. 


(3.10) 


If  we  specialize  to  the  case  (2.25),  and  assume  that  the  market  price  Vt  follows  a  mean 
reverting  process  with  constant  coefficients,  then 


■^'(^T  I  Po)  = /(?7o»voiT)  (3.11) 

where  go  is  observed  and  Vq  still  needs  to  be  inferred  from  traded  options.  This  is  only  a  two- 
dimensional  quantity,  however,  so  only  two  expiration  dates  are  needed.  Note  that  T  now  represents 
the  observed  filtration. 

For  an  explicit  expression  for  /,  suppose  that 


Vt  =  VL  +  Vu 


(3.12) 
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where  vi  is  the  long  run  value  of  Also  assume  that  Vt  follows  a  mean  reverting  process 


dvf  —  Aifdt  “h  AodZt, 


(3.13) 


Then 


and  so 


f{no,vo,t)  = 

y  TM~^^'^Vsds  +  j  I  ^o)  ,  (3.14) 


f{f]Q,  Vo,  t) 


=  (l,0)(t/') 


—  1 

e"  '7?o 


+  (1,0)  [(t/')“^(e'''‘  -  u't  -  vl 

+  (1,0)  [Ut  +  (Af^/2  _  A-\e^^  -  I) 


(vo  -  vl) 


where  is  the  solution  of 


(3.15) 


i^'Qt  -  ^tA  =z  (i/')“'rM-'/2e^'  +  (t/')-'e^''r.¥-i/2  ^ 

This  joUy  little  expression  is,  of  course,  more  suitable  for  software  than  for  human  consumption, 
but  the  point  is  that  all  the  quantities  in  (3.15),  except  for  fjo  and  uq,  are  practically  observable 
in  the  long  run.  fjo  is  observed,  and  one  only  needs,  therefore,  two  implied  quantities,  namely  the 
two  components  of  uq.  This  wiU  use  two  expiration  dates  to  infer  the  relevant  quantities.  If  there 
are  more  dates  available,  one  can  use  regression,  and  one  can  also  do  model  checking. 

How  is  fjQ  computed?  For  the  purpose  of  our  limiting  argument,  any  predictor  fj^  which 
converges  in  probability  to  fjo  as  6  tends  to  zero  will  do.  Further  fine  tuning,  as  in  Florens-Zmirou 
(1993),  Chesney,  Elliot,  Madan  and  Yang  (1993)  and  Pastorello  (1996),  is  obviously  desirable  in 
practice,  though  we  shall  not  pursue  this  issue  here. 
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4.  CONCLUSION 

We  have  seen  in  the  above  that  asymptotic  methods  can  be  effective  in  studying  incomplete 
markets.  Also,  it  is  apparent  that  it  is  undesirable  to  rely  on  implied  quantities  only,  and  that  one 
needs  to  also  use  historical  data  to  extract  information  about  unobserved  quantities.  The  effect  of 
estimation  is  incorporated  into  the  options  price  in  a  natural  manner. 

The  specific  model  and  the  specific  type  of  asymptotics  used  is  by  no  means  meant  to  be 
anything  more  than  a  test  case  and  an  example.  There  is  room  for  substantial  exploration  into 
how  one  can  best  implement  the  methodology  advocated  here. 
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APPENDIX 


Proof  of  Theorem  1.  Write 


log  sf  =  log  5o  +  /  Hudu  +  if 
Jo 

[  {<tI  -  al)du  +  L\ 

Jo 

and 

ni  =  I\\u  -  6-^\L)du  +  L^. 

Jo 

L  =  {L^,L^,  L^)  is  then  a  martingale  with  predictable  covariation  matrix  given  by 

r  ] 

{L,L)t= 

{L^,L^)t  Sf^Xads, 

'  (T\t  k^t  0 

=  kst  kit  alt  4-Op(l) 

0  cr?t  A/,t_ 


(A.l) 


Since  Esup(A£f)‘*  =  o(^),  and  since  L  and  W  have  zero  covariation,  it  follows  from  Theo¬ 
rem  VIII.3.12  (p.  432-433)  in  Jacod  and  Shiryaev  (1987)  that  {L,  W)  converges  in  law  to  a  Brownian 
motion,  with  L  independent  of  W  and  with  covariance  matrix  given  as  the  limit  in  (A.l).  The 
result  then  follows  since  /q  —  <^i)ds  and  ^^/^(As  —  6~^\[^)ds  converge  in  probability  to 

the  integrals  of  the  respective  limits,  by  the  uniformity  in  (2.2)-(2.3).  Q.E.D. 

Proof  of  Theorem  2.  By  Ito’s  Lemma  (see,  e.gr.,  Theorem  1.4.57  (p.  57)  of  Jacod  and 
Shiryaev  (1987)), 


=  AC-{S,‘, I)  -  <)A5; 

-  \C-ss(Sl,t)(Slfaldt 

A 
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where  we  have  used  that  C*  +  rC^S  +  \C*sgS^cT^  =  0.  The  dependence  on  is  suppressed  in  the 
notation.  By  Taylor  expansion, 

at  =  i  f  Ch{SL,«)d(i 

Jo 

+  f{Zlu)d[logS^logS\logS%  (A.3) 

where  |Zf  —  5fl  <  |A5f|.  Z  is  measurable  and  cadlag  by  construction,  and  by  Lemma  VI.3.31 
(p.  316)  of  Jacod  and  Shiryaev  (1987)  and  (our)  Theorem  1,  Z^  converges  to  S  in  law.  By  the  con¬ 
tinuity  assumption  of  the  current  theorem,  f{Zf,-)  converges  to  in  law  (Remark  VI.3.8 

(p.  312)  in  Jacod  and  Shiryaev  (1987)).  One  can  then  use  Theorem  VI.3.21  (i)  (p.  314)  in 
the  same  work  together  with  uniform  integrability  to  conclude  that  the  second  term  on  the 
r.h.s.  of  (A.3)  converges  as  specified  to  the  second  term  on  the  r.h.s.  of  (2.16).  This  is  because 
i“^/^[log5^,log5^,log5'^](  =  [^^,log5'^]j,  the  convergence  properties  of  which  follow  from  Corol¬ 
lary  VL6.7  (p.  342)  in  Jacod  and  Shiryaev,  in  view  of  (our)  Theorem  1. 

The  first  term  on  the  r.h.s.  of  (A.3)  is  handled  in  analogy  to  the  proof  of  Theorem  1.  Q.E.D. 

Proof  of  Theorem  3 — Solution  of  Equation  (3.4).  Since  the  solution  is  obvious  if 
u'  is  symmetric,  we  shall  exclude  this  case  in  the  following. 

Set  r  =  u'  —  and  7  =  M^/^7.  The  problem  is  equivalent, 

obviously,  to  finding  a  unique  symmetric  nonnegative  definite  (u.s.n.d.)  solution  to 

(F  -  ?')(f  -  u'y  =  v'{u'y  -1-  77*.  (A.4) 

If  one  writes  the  right  hand  side  of  (A.4)  as  V LV*  where  L  is  diagonal  and  V  is  orthonormal,  one 
can  further  set  f  =  V*TV  and  D'  =  V^u'V ,  and  the  problem  then  reduces  to  finding  a  u.s.n.d. 
solution  to 

(f  -  P'){f  -  D'y  =  L.  (A.5) 

The  set  of  T  satisfying  (A.5)  with  no  further  constraint  is  given  by 

(A.6) 


f  =  -I-  9' 
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where  U  can  be  any  orthonormal  matrix. 

Since  f  must  be  symmetric,  the  two  off  diagonal  elements  of  +  u'  must  be  the  same, 


=  I'u  -  ^21- 


Also,  since  U  is  orthogonal, 


U11W21  +  U22«12  =  0. 


Solving  these  two  equations  yield  that 


where 


U2I  =  U22{l^u  -  ^'2i)/det 

U12  =  -un(i>i2  -  i^2i)/det. 


det  = 


(we  discuss  below  why  det  ^  0).  +  tij2  =  1  is  equivalent  to 


and,  similarly,  “21  +  y>22  ~  1  to 


(A.8) 


(A.IO) 


(A.11) 


(A.12) 


Hence  un  =  ±xi22*  The  minus  sign  is  ruled  out.  however,  since  otherwise  (by  (A. 9))  U  would  be 
symmetric  and  hence  (since  orthonormal)  the  identity.  This  could  mean  that  P[2  =  ^211  whence  i/' 
would  be  symmetric,  which  we  have  previously  ruled  out. 

A  solution  exists,  therefore,  provided  (A. 11)  and  det  7^  0  (if  det  =  0,  there  is  no  solution, 
since  the  case  u'  symmetric  has  been  ruled  out).  This  is  equivalent  to 


(A.13) 


This  condition  is  satisfied  since 


(1/(2  -  1^21  <  tr(i/V'*) 
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<  tr(i) 

(A.14) 


There  are  now  two  possible  solutions  of  our  system,  corresponding  to  positive  or  negative 
sign  of  uii.  However, 

tr(f )  =  uu{lT  +  4^')  +  (A.15) 

which  is  negative  if  un  <  0.  This  contradicts  the  requirement  that  f  be  positive  definite.  Hence, 
the  only  possible  solution  for  U  is 


U 


\  a  (1  — 


where 


i?i2  -  i>' 


a  = 
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(A.16) 


(A.li 


This  solution  exists  (z.e.,  \a\  <  1),  T  satisfies  (3.4),  and  is  symmetric.  It  is  also  the  only  solution 
which  may  be  positive  definite.  To  show  that  this  requirement  is  indeed  satisfied,  with  F  =  QAQ*. 
and  note  that  A  can  be  chosen  so  that  F  is  positive  definite,  and  hence  this  must  be  a  property  of 
the  solution  we  have  described.  Q.E,D. 


Proof  of  the  Remainder  of  Theorem  3.  Provided  (i)  holds,  it  follows  from  Theo¬ 
rem  12.7  (p.  33)  of  Lipster  and  Shiryaev  (1978)  that  F  is  the  solution  of  (3.4).  Equation  (3.5)  is  a 
consequence  of  the  same  result. 

(i)  ^  (ii)  follows  from  the  above  and  from  Theorem  12.6  (p.  31)  of  Lipster  and  Shiryaev 
(1978). 

(ii)  (i);  Let  P  be  the  probability  distribution,  and  let  Q  be  one  in  (i).  By  our  assumptions. 
dP/dQ  exists  and  is  J^Q-measurable.  Since  i/'  is  negative  definite,  (rjt)  is  ergodic,  and  so  rjt  is 
asymptotically  independent  of  P'q  and  hence  dP/dQ»  It  follows  that,  for  bounded  continuous 


E<}  (</(?,)§ 

Eq  1 

(  dP 
[dQ  ' 

^obs^ 

1 
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=  £e(s(>;.)l^“)+Mi)- 


Hence,  by  (ii),  we  must  have  Q  =  P,  and  (i)  follows. 


Q.E.D. 
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